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ABSTRACT 

Self-interacting dark matter has been proposed as a hypothesis to explain the shallow central slopes of 
the density profiles of dark matter halos in galaxies. In order to be consistent with observational studies 
at scales of galaxy clusters, the cross section should scale inversely with the velocity of collision. In this 
paper we consider the mass density profile of the halo of the low surface brightness (LSB) galaxy NGC 
5963 to place an upper limit on the dark matter cross section for collisions with velocities ~ 150 km s~\ 
i.e. at the low velocity regime. After calibrating against cosmological simulations, we found that the 
large inferred dark matter concentration and central dark matter density in NGC 5963 are inconsistent 
with an effective collisional cross section per unit of mass > 0.2 cm 2 g -1 . Corrections were applied 
in order to account for reduction of the core by the adiabatic contraction caused by cooling baryons. 
Our limits that involve a number of simplifying, but always conservative, assumptions, exclude the last 
permitted interval for velocity-dependent cross sections to explain the flat density core in LSB galaxies. 
Implications for the nature of dark matter are also discussed. 

Subject headings: dark matter — galaxies: halos — galaxies: kinematics and dynamics — galaxies: 
spiral 



1. INTRODUCTION 

Numerical studies of structure formation with the colli- 
sionless cold dark matter (CDM) scenario predict dark ha- 
los with steep central cusps (e.g., Navarro, Frenk & White 
1996, 1997, hereafter NFW), whereas most of the rota- 
tion curves of dwarf galaxies and low surface brightness 
(LSB) galaxies suggest that their halos have constant den- 
sity cores (e.g., Marchesini et al. 2002; de Blok & Bosma 
2002, and references therein). Self- interacting dark matter 
(SIDM) with cross sections per unit of mass in the range 
0.5-6 cm 2 g _1 was proposed as a possible route to reduce 
the density cusp, as thermal conduction replaces the cen- 
tral cusp by a soft core (Spergel & Steinhardt 2000). A 
tight constraint on the cross section just under 0.1 cm 2 
g _1 has been derived by considering the formation of gi- 
ant cluster arcs (Meneghetti et al. 2001) and the size of the 
cores of galaxy clusters (e.g., Arabadjis, Bautz & Garmire 
2002; Lewis, Buote k Stocke 2003; Arabadjis & Bautz 
2004). Since this cross section is too small to produce 
galactic cores, it has been pointed out that a velocity 
dependent cross section might reduce the effects of self- 
interaction on cluster scales (e.g., Firmani et al. 2000). 
For a cross section per unit of mass, which varies as some 
power of the relative velocity between colliding particles 
o"dm = c*(wo/w r oi) a , where vq = (100 kms -1 ), there exists 
a certain range of values of er* and a for which the halos of 
dwarf galaxies should present a core, while halos of galaxy 
clusters would not change their core sizes, ellipticities and 
arcs significantly. 

Constraints on the parameters cr* and a were obtained 
by requiring that dwarf galaxies observed today have yet to 
undergo core collapse and that dark halos must survive the 
heating from hot cluster halos (Gnedin & Ostriker 2001; 
Hennawi & Ostriker 2002). All these requirements are 
satisfied for a very narrow range of parameters cr» = 0.5- 



1 cm 2 g _1 and a = 0.5-1. For cross sections within this 
suitable range, flat cores with densities of ~ 0.02 M Q pc -3 
are formed in the central regions of galactic halos as it was 
confirmed numerically by Colin et al. (2002) in cosmolog- 
ical simulations of SIDM. 

Low surface brightness galaxies are ideal to put upper 
limits on the strength of dark matter (DM) self-interaction 
at the low- velocity regime (i.e. relative velocities of ~ 150 
km s^ 1 ). In contrast to relatively recently formed objects, 
like clusters of galaxies, that present unrelaxed mass dis- 
tributions, LSB galaxies with high central densities may 
have the highest formation redshifts and, therefore, they 
have had almost a Hubble time to soften their cusps. In 
the present study we concentrate on the implications of 
the halo of the LSB galaxy NGC 5963 for dark matter 
self-interaction. The high central density and the small 
core radius of the halo of NGC 5963 suggest that either 
the DM cross section, at the low velocity regime, is rather 
small < 0.1 cm 2 g _1 or the dark halo is undergoing an un- 
desirable, dramatic core collapse. Therefore, it seems very 
unlikely that collisional scattering between DM particles 
is the main agent for the formation of cores. 

We will start in §2 with a description of the basic physics 
behind the scenario of SIDM and some predictions, which 
will be used in the remainder of the paper. In §3 we briefly 
describe the properties of NGC 5963 and present mass 
models over a range of mass-to-light ratios of the stellar 
disk. We then compare with predictions of cosmological 
SIDM simulations to constrain the self-interaction cross 
section of DM particles (§4). Some implications of the 
results are discussed in §5. 

2. SELF-INTERACTING DARK MATTER AND THE SIZE OF 
THE CORES 

The mass-density profiles of numerically simulated colli- 
sionless CDM halos are commonly parameterized with the 
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analytical Navarro, Frenk & White profile: 



where r s is the characteristic radius of the halo and pi is 
related to the density of the Universe at the time of col- 
lapse. Within the core region, the dark matter density 
increases as a power law p <~ and the velocity dis- 
persion decreases toward the center at the radius r s . The 
parameters {pi,r g ) are usually expressed in a slightly dif- 
ferent form by the concentration parameter c = r2oo/?"s 
and the rotation velocity V200 at radius r 2 oo- The latter 
is the radius inside which the average overdensity is 200 
times the critical density of the Universe. 

If dark matter particles experience self-interaction, the 
scattering thermalizes the inner regions of dark halos, pro- 
ducing a constant-density core (e.g., Burkert 2000). It is 
expected that the core radius will be comparable to the 
inversion radius, i.e. the radius at which the velocity dis- 
persion peaks. After reaching the maximum core size, the 
core can begin to recollapse due to the gravothermal catas- 
trophe. In order to have a core, the relaxation time should 
be less than <~ 1/10 of Hubble time, so that it should have 
sufficient time to form a core by thermalization, while the 
lifetime of its core should exceed 10 Hubble times, so the 
core should be far from collapse. 

In the monolithic scenario in which the halos evolve at 
isolation, the core radius increases in time as the result 
of self- interaction between halo particles. Starting with a 
profile with inversion radius , the core radius of the pseu- 
doisothermal profile, r c , achieves a radius w 0.4rj in one 
core-radius relaxation timescale t rc = (3/9CTd m w) _1 , with 
p the density and v the velocity dispersion, both evalu- 
ated at 0.6r; (e.g., Kochanek & White 2000) 1 . The maxi- 
mum radius that the core reaches before recollapse may be 
even larger. In fact, for the particular case of a Hcrnquist 
(1990) model with break radius m, Kochanek & White 
(2000) found that r c reaches a maximum value of 0.23ry. 
Since the velocity dispersion for a Hernquist profile peaks 
at 0.33rf/, the core radius r c at the maximum expansion 
is approximately 0.7^. 

The assumptions behind the monolithic scenario that 
the initial configuration follows a NFW profile and that 
halos evolve at isolation appear to be inconsistent with the 
SIDM hypothesis in a cosmological context. Cosmological 
simulations of the formation of SIDM halos, either with a 
constant cross section or dependent on velocity, have been 
carried out by Yoshida et al. (2000), Dave et al. (2001), 
Colm et al. (2002) and D'Onghia, Firmani & Chincarini 
(2003). The inclusion of dynamically hot material that 
continuously is accreted onto the halo (in a cosmological 
context) can prevent efficient heat transfer from the core 
to the halo. Moreover, it turned out that since collisions 
already modify halo profiles in the dynamical process of 
halo growth, there is a smooth trend of increasing core ra- 
dius with the effective value of <7dm- A core with a radius 
<~ 0.4r s is formed if the mean number of collisions per par- 
ticle and per Hubble time in the center halo is about 4-6. 
We will use this result in §4. 

The existence of galaxies with large DM concentrations 
and central DM densities implies upper limits on the num- 

1 Kochanek & White (2000) estimated the core radius by the point at 
have used that r c is related to r 1 / 4 by the relation r c = r 1 / 4 /l.7. 



ber of collisions per particle and hence on a„. In the fol- 
lowing we consider the implications of the concentrated 
halo in the galaxy NGC 5963. 

3. THE GALAXY NGC 5963 AND ITS DARK HALO 

3.1. NGC 5963 and its Rotation Curve 

NGC 5963 is a relatively isolated LSB galaxy associated 
with the NGC 5866 group, with a suitable inclination an- 
gle of 45 ± 4 deg. Its adopted distance is ~ 13 Mpc, based 
on the recession velocity derived by Bosma, van der Hulst 
& Athanassoula (1988) and a Hubble constant of 70 km 
s _1 Mpc -1 . The uncertainty in the distance is probably of 
order 3 Mpc. The distributions of light and mass were first 
reported by Romanishin et al. (1982), who found a lumi- 
nosity profile that exhibits a lens component, i.e. a com- 
ponent with a plateau and a steep outer edge, and a small 
bulge (see also Simon et al. 2004) . The extent of the lens is 
about 1.4 kpc. The Hi rotation curve was constructed by 
Bosma et al. (1988). These authors emphasized that the 
dark halo of this galaxy is more concentrated than in nor- 
mal Sc's with similar rotation curves. In this sense, NGC 
5963 obeys the rule pointed out by Sancisi (2003) that 
an excess of light corresponds an excess of rotation and 
viceversa. Beyond a radius of 7 kpc the rotation curve 
becomes fairly flat with a circular velocity v c ~ 131 km 
s _1 (see Fig. 1). If the dark halo can be approximated by 
the pseudoisothermal sphere, the one-dimensional veloc- 
ity dispersion for halo particles for NGC 5963 would be of 
v s=s v c /V2 = 90 km s _1 and, therefore, the halo particles 
undergo collisions in the low-velocity regime. 

In order to study the cuspiness of the mass density pro- 
file in the galaxy centers it is required a rotation curve with 
a high spatial and velocity resolution and a large exten- 
sion. This can be achieved by combining CO, Ha data, and 
radio H 1 data. NGC 5963 is one of the targets of the sam- 
ple of galaxies with high-resolution two-dimensional Ha 
and CO velocity fields, which are free of beam-smearing, 
in Simon et al. (2004). Two-dimensional velocity fields 
are useful in order to avoid systematic uncertainties and 
to determine the existence of radial motions. Multiple 
wavelengths studies of the velocity fields may also help 
to elucidate the origin of systematic errors that can make 
density profiles appear artificially shallow. The observa- 
tions in CO and Ha provide the rotation curve within a 
galactocentric radius of 4 kpc, whereas the Hi rotation 
curve allows us to trace the potential out to 10 kpc. This 
galaxy lacks of measurable non-circular motions, i.e. ra- 
dial velocities < 5 km s _1 , except perhaps in the rings 
between 12" to 30", those being not large enough to af- 
fect the derivation of the DM density profile. Simon et 
al. (2004) compared the Ha and CO velocities at every 
point and found an excellent agreement between them, 
with a mean offset of less than 1 km s _1 , and a scatter 
of 7.8 km s _1 . Though the uncertainties in the rotation 
curve can be reduced by combining both velocity data, we 
prefer to use the CO rotation curve and its associated er- 
rorbars rather than the formal errorbars of the Ha rotation 
curves, which simply indicate that the Gaussian fit to the 
profile was well-determined. In particular, the CO rota- 
tion curve is less sensitive to asymmetric drift corrections 

which the density drops to 1/4 the central density. To estimate r c , we 
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and other magnetic terms, due to the smaller small-scale 
motions of the molecular clouds compared to the ionized 
gas (Sanchez-Salcedo & Reyes- Ruiz 2004). In addition, 
other uncertainties in mass models, as the mass-to-light 
ratio of the different components, impose a minimum er- 
ror in any fit. Figure 1 shows the CO rotation curve in 
the inner galaxy and the Hi rotation curve in the outer 
parts, with error bars indicating the uncertainties in the 
observations. 

The present galaxy is an ideal laboratory with which to 
test the effects of collisions between DM halo particles in 
the low velocity regime, not only for the properties of the 
galaxy but also for the quality of the data. 



3.2. Mass Models and Adiabatic Contraction 

We have constructed different mass models for this 
galaxy. First we calculate the rotation curve due to the 
neutral (H i+Hc) gas mass from the radial distribution of 
the Hi surface density from Bosma et al. (1988). The sur- 
face brightness profiles in the optical and near-/i? bands 
given by Simon et al. (2004) were used to derive the disk 
mass distribution by adopting a constant mass-to-light ra- 
tio for the disk. The resulting mass distribution was con- 
verted into a disk rotation curve. The difference between 
the sum of the squares of the observed velocity and the 
stellar plus gaseous rotation velocities was converted into 
a halo mass distribution after assuming a certain form for 
the halo. Along this paper, all the fits to the CO rotation 
curve were performed by a ^-minimization. To make the 
best use of both the CO and H I data, we have used a hy- 
brid rotation curve. It consists of the CO data over the 
range of radii where available and 21 cm data to define 
the outermost points. This was accomplished by minimiz- 
ing x 2 between and 4 kpc but requiring that between 5 
and 10 kpc, the goodness of the fit must be better than 
la level, while rejected otherwise. In this way, the fits are 
statistically consistent with the H I rotation curve at large 
radii. 

Different mass decompositions were fitted for the stellar 
mass-to-light ratio in the i?-band varying from (mini- 
mum disk) up to 1.2 Mq/Lq (maximum disk). We will de- 
note by the stellar mass-to- light ratio in the i?-band in 
solar units. It is noteworthy that T* = 1.2 is the maximum 
value allowed by the smooth rotation curve but avoiding 
a hollow halo in the inner two points. However, it is sub- 
maximal in the sense that the maximum contribution of 
the stellar disk to the rotation curve is less than 0.8 times 
the observed maximum rotation. In practice, the stellar 
curve is scaled until the inner points match those of the 
smooth curve. Depending on how many points are used, 
the maximum disk mass-to-light ratios may vary as much 
as <~ 15%. In the case of NGC 5963 and using the inner 
five points, Simon et al. (2004) obtained a dynamical value 
for the maximum mass-to-light ratio larger than the value 
quoted above, but it appeared unrealistic and inconsistent 
with the predictions from the galaxy colours by the Bell 
et al. (2003) population synthesis models. The dynami- 
cal maximum disk inferred by using two points is fully 
consistent with these predictions. 

Fig. 2 shows the best-fitting mass models for different 
values of Y* with both the NFW halo and the cored pseu- 



doisothermal halo: 

P'W = -^TT2- ( 2 ) 

Primed quantities are used where appropriate to denote 
the (observed) final states after adiabatic compression of 
the dark halo by the infall of the baryons as they cool 
and settle into a disk. The parameters of each model, the 
reduced x 2 of the fit and the probability p that the data 
and the model could result from the same parent distribu- 
tion are given in Table 1. It is important to remark that 
the observation that pseudoisothermal fits produce low re- 
duced x 2 values does not demonstrate the presence of a 
constant-density core. Simon et al. (2004) have shown 
that a power-law density profile p cx r -12 fits also very 
well the rotation curve of NGC 5963. 

In order to compare with simulations which do not in- 
clude the baryonic physics, we need an estimate of the 
DM profile before the adiabatic compression by the baryon 
condensation, i.e. the unprimed parameters of the halo. 
A cored profile is expected after adiabatic decompression. 
Since the core baryons cause the dark profile to steepen, 
we have fitted the observed rotation curve, for a given 
Y* > 0.2, with a cuspy profile (for convenience we have 
chosen the NFW profile) and recovered the halo parame- 
ters (po> r c) of the approximately pseudoisothermal sphere 
before compression, following the standard model (e.g., 
Blumenthal et al. 1986; Flores et al. 1993). For Y* < 0.2, 
the effect of compression by the baryons is so small that 
a cored profile is assumed before and after contraction. 
The recipe of Blumenthal et al. (1986) assumes that the 
matter distribution is spherically symmetric and that par- 
ticles move on circular orbits. Jesseit et al. (2002) veri- 
fied that if the disk formation is smooth, the contracted 
halos have density profiles that arc in excellent agree- 
ment with the standard predictions. If the disk forma- 
tion involves clumpy, cold streams or if the bulk of the 
central stars formed when the halo mass was still being 
assembled, the level of halo contraction may be smaller. 
Gnedin et al. (2004) suggested an improvement of the stan- 
dard model of the response of the halo to condensation of 
baryons by a simple modification of the assumed invari- 
ant (see Gnedin et al. 2004 for details). They show that 
the standard model systematically overpredicts the con- 
traction. Moreover, there exist other empirical evidence 
which suggests that adiabatic contraction is avoided (e.g., 
Loeb & Peebles 2003; Dutton et al. 2005). Therefore, our 
procedure provides an upper limit on the effect of adiabatic 
contraction. 

For illustration, Fig. 3 shows the contribution to the ro- 
tation curve of the dark halo after correcting by adiabatic 
contraction, in the most extreme case of maximum disk 
Y* = 1.2. The fit to it adopting the pseudoisothermal 
profile is also shown. Note that if we demand that the 
density of the halo must decrease monotonically with in- 
creasing radius, values Y* > 1.2 would not be permitted 
in the standard adiabatic decompression approach. 

In Table 1, best-fitting halo parameters are reported as a 
function of T*. Even adopting the maximum disk solution, 
the halo of NGC 5963, at present time, is extraordinary 
concentrated, with a central dark matter density greater 
than 0.35 M Q pc~ 3 and a core radius < 1 kpc. A thorough 
discussion aimed to identify what about NGC 5963 makes 
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it unique was given by Simon et al. (2004) and need not be 
repeated here. In the next sections, we will discuss the im- 
plications of the halo of NGC 5963 for the self-interaction 
of DM particles. 

4. CONSTRAINTS ON THE CROSS-SECTION OF DM 
PARTICLES AT THE LOW VELOCITY REGIME 

In principle we do not really know whether the core of 
the galaxy under consideration is still in the expansion 
phase or undergoing gravothermal core collapse aided by 
the gravitational contraction due to the baryons. In fact, 
the core contraction of the halo by the collapsing baryons 
may trigger a fast core collapse (e.g., Kochanek & White 
2000, and §4.3). This effect may be important for T* ra- 
tios close to the maximum disk value. Using cosmological 
simulations of SIDM as a calibrator, we first put limits 
on the cross section of DM particles assuming that the 
core of NGC 5963 is in the expanding phase due to heat 
transfers inwards (§4.1). For large values of T*, however, 
core collapse triggered by baryon condensation may have 
dramatic consequences for the longevity of the core. We 
have to demand that the halos are in no danger of collaps- 
ing. This requirement implies another bound on the cross 
section, which is discussed in §4.3. 

The galaxy halo profiles and the tight cross section lim- 
its coming from clusters of galaxies can be reconciled only 
if the cross section were inversely proportional to the halo 
velocity dispersion. Nevertheless, except for subhalos un- 
dergoing the action of the hot halo environment, most of 
the relevant results found in simulations with a cross sec- 
tion independent of halo velocity dispersion are valid if the 
value of the cross section is interpreted as that appropri- 
ate to the velocity dispersion of the halo being considered. 
Therefore, since we are primarily interested in a single ob- 
ject, we can gain additional physical insight by making 
also use of simulations with a velocity-independent cross 
section. 

4.1. Halos in Midsized Galaxies under SIDM: 
Calibrating Collisional Effects 

In cosmological simulations of structure formation in 
SIDM scenario with effective cross sections 0.5-1 cm 2 g _1 , 
midsized halos (masses <~ 10 10 ~ n M©) present cores with 
radii in the range r c w 2.5-5 kpc and central densities 
po w 0.01-0.06 M© pc" 3 (Yoshida et al. 2000; Dave et 
al. 2001; Colin et al. 2002; D'Onghia et al. 2003). For 
instance, D'Onghia et al. (2003) derive r c = 5-6 kpc and 
Po ~ 0.02 M© pc~ 3 for a halo with maximum circular ve- 
locity of 120 km s _1 . Therefore, effective cross sections 
in the range 0.5-1 cm 2 g _1 are suitable to reproduce soft 
cores in late- type galaxies 2 . 

All the simulations with a^m = &*(vo/v) also shown 
that for cr* = 0.5-1 cm 2 g -1 , the number of collisions per 
particle per Hubble time at the center of the halos, was be- 
tween 4-6, roughly independent of the halo mass. In fact, 
while a few collisions (between 2 and 3) at the halo center 
are enough to produce a constant density core (Yoshida 
et al. 2000), 4-6 collisions are sufficient to produce a core 
with r c ~ 0.4r s0 , with r s0 the scale radius of the NFW 



halo when it is rcsimulated in the standard, collisionless 
case ((7dm = 0). The reason is that each scattering pro- 
duces a change Av ~ v, and particles escape the core in a 
single scattering. 

Rescaling the number of collisions per particle at the 
center of the cluster labelled SlWa in Yoshida et al. (2000), 
we find 



N col 



Po 



0.02 M Q pc- 3 J VlOOkms 



1 cm 2 g 



(3) 

We stress here that po is the final unprimed central den- 
sity of the core, after thcrmalization. Let us estimate AT co i 
in the halo center of NGC 5963 for mass models with 
T* < 0.7. For a DM cross section <7d m '5ioo between 0.5 and 
1.0 cm 2 g _1 , where v W q = (u/100kms -1 ), and a central 
density po w 0.4 M© pc -3 , the mean number of collisions 
is iV co i ;=y 20-45. Hence, we expect NGC 5963 having a 
core radius greater than 0.4r s o- In fact, from the scaling 
laws of Yoshida et al. (2000), we see that the core radius 
reaches a size ~ 1.5r s o when the collision rate per parti- 
cle at halo center is ~ 30 per Hubble time. In order to 
quantify the expected size of the core of this galaxy under 
SIDM, we need its characteristic radius r S Q . 

An estimate of r s o can be inferred for NGC 5963 by fit- 
ting the rotation curve beyond 5 kpc with a NFW profile. 
The observed inner rotation curve cannot be included as 
it is altered as a consequence of DM self-interactions and 
the adiabatic contraction that causes the DM profile to 
steepen. Since there are many combinations of the pa- 
rameters (c, V200) that can fit the rotation curve, we need 
to choose either c or V200 • Expected values of the concen- 
tration in ACDM cosmology have a 2a range from 5 to 
22 (Eke et al. 2001; Jing & Suto 2002). We will take a 
value c = 20, which lies among the largest values in the 2a 
uncertainty. The exact value of r s o also depends slightly 
on T*. For the intermediate case T* = 0.6 (and c = 20), 
we find that the characteristic radius is r s o ~ 6 kpc. For 
the average value predicted by simulations, c = 13, we 
found r s0 = 11.7 kpc. Hence, as a conservative calibra- 
tion, we will assume that for A co i > 6, a core radius 
r c :=y 0.4r s o > 2.5 kpc should have been formed. 

In Table 1 we see that for models with T* < 0.7, the 
inferred core radius before contraction in NGC 5963 is 
smaller than 1.2 kpc; the question that arises is how many 
collisions suffice to develop a core with r c < 1 kpc. To an- 
swer this question we have examined the SIDM simulations 
of the evolution of a Hcrnquist halo with total mass Mt 
and break radius rn in Kochanek & White (2000). These 
authors carried out simulations surveying the dimension- 
less cross section a<\m = -MrCdmAff • According to their 
Fig. 2c, the core density in the phase of core expansion, and 
for halos with different <7d m , goes as cx er^, with r\ < 1/2. 
Moreover, at the time that the core of the simulation with 
dimensionlcss cross section <7d m = 3 reaches its maximum 
radius, the run <7d m = 0.3 has already developed a core of 
half this radius, with a mean collision count at halo center 
<~ 5 times smaller. Applying this scaling to NGC 5963 it 
holds that a mean of 6/5 collisions per particle at the cen- 
ter, will suffice to produce a core of 2.5/2 kpc. Although 

2 In order to compare simulations, we notice that if the cross section is assumed to be inversely proportional to the collision velocity, a^m = 
cr* (uo /v Te \ ) , the effective cross section for a halo with (one-dimensional) velocity dispersion v is <t*/(v / ?™ioo), where vioo = (C/100 km s _1 ). 
To show this relation one has to remind that, for a Maxwellian velocity distribution, (l/v re i) = l/(y^rrv). 



5 



these estimates are based on simulations of the relaxation 
of an isolated galaxy having initially a cuspy Hernquist 
profile, it was also found in cosmological simulations that 
a few collisions (exceeding 2) suffice to develop a kpc-size 
core (e.g., Yoshida et al. 2000). To keeps matter simple 
we will adopt the generous condition that only if iV co i < 2 
then r c < 1 kpc. 

Putting together, we take the following, rather conserva- 
tive, relationships as a calibrator of the collisional effects: 



if N col = 6 <=> r c = 2.5 kpc, 
if 7V co i = 2 r c < 1.0 kpc. 



(4) 
(5) 



Suppose that (7d m cx v^ 1 . The generalization for a power- 
law <7dm v~ a is obvious and it will be ignored in the 
interest of simplicity. To place an upper limit on CTdm^ioO) 
we proceed as follows. We first derive the best-value fit- 
ting parameters (r c , po) for a given T*. If 1 kpc < r c < 3 
kpc, we estimate N co \ by interpolating Eqs. (4)- (5) linearly. 
Once iV C oi is known, Equation (3) immediately provides us 
with an upper limit on cr = Cdm^ioo ■ ^ r c i s smaller than 
1 kpc, we will assume that N co i < 2, which gives: 

*o < OJcmV 1 ' P ° 



^O.O2M pc- 3 y ' ^ 

For instance, if we have a situation in which r c < 1 kpc 
and po — 0.4 M Q pc~ 3 , we would obtain er < 0.035 cm 2 

4.2. Results 



Figure 4 shows the resulting upper limits on Co versus 
Y*, following the procedure described in the previous sub- 
section. Confidence levels were calculated as a measure of 
the sensitivity of the results to the goodness of the fit to 
the rotation curve. Other sources of uncertainties in the 
parameters of the model, as the distance to the galaxy, its 
inclination, or the adopted characteristic radius r s n, were 
not taken into account in these contours. We must warn 
that these error bands cannot be considered as real proba- 
bility indicators because the velocities and their errors are 
not free of systematic effects. 

Our derivation overestimates the upper limit on cr be- 
cause we have assumed in the derivation of the halo pa- 
rameters that baryons only act to hasten core contraction. 
Nevertheless, other phenomena as galactic bars, outflows, 
massive black holes in the halos, and other processes as- 
sociated with AGN or star formation have been proposed 
as mechanisms to erase cuspy halos. We are being conser- 
vative in ignoring these phenomena because incorporating 
their effects would only serve to place more stringent con- 
straints on (To. 

The bound on cr depends strongly on Y*. From Fig. 4 
we see that models with Y* w 0.6 are consistent with <tq 
just under 0.1 cm 2 g _1 . For T* = 0.35, a > 0.05 cm 2 
g _1 can be excluded virtually at 2a confidence interval. 
Within the range < T* < 0.7, our limit excludes the in- 
terval proposed to explain the flat mass profiles in galaxies. 

For Y* > 1.0, the present analysis only disfavours the 
interval of astrophsical interest, cr = 0.3-1 cm 2 g _1 , at 
less than 68% confidence level. We should notice, how- 
ever, that if the contribution of the baryons to the po- 
tential well is important, the core may be in danger of 

3 It is simple to show that o^m and (A//f) e( j ge are related by the relation 



undergoing gravothermal collapse. The requirement that 
the core is not undergoing a fast collapse will provide a 
more stringent upper limit for the DM cross section for 
T* > 0.8. 



4.3. Core Collapse 

In the latter section we found that a cross section 
(To ~ 0.5 cm 2 g -1 and Y+ > 1.0 would be marginally con- 
sistent with a picture in which the halo of NGC 5963 might 
have formed a core with a radius ~ 2-2.5 kpc, although 
significantly reduced because of the adiabatic contraction 
by baryons. Nevertheless, we must also require that this 
galaxy has yet to undergo core collapse. For these high Y*, 
the cooling baryons will compress the dark matter and ac- 
celerate the core collapse in different ways (e.g., Kochanck 
& White 2000). By raising the density due to the contrac- 
tion, the relaxation time drops. In addition, the adiabatic 
compression produces a steeper central density cusp which 
further speeds up the ultimate evolution of the system. 
Moreover, the initial inversion of the temperature (veloc- 
ity dispersion of DM particles) profile responsible for the 
core expansion may be erased by the adiabatic heating of 
DM particles, leading directly to a core collapse with no 
room for the expansion phase. We need to estimate the 
characteristic time for core collapse in the case Y* > 0.8, 
i.e. when the effects of contraction are important. In their 
numerical simulations, Kochanek & White (2000) assumed 
as initial conditions a Hernquist profile with a cuspy core. 
In our case, i.e. for values Y* > 0.8, it is realistic to start 
with such an initial configuration because a cusp in the 
density profile is a natural consequence of the adiabatic 
contraction caused by the baryons. In this regard, we can 
make use of the quantitative study of Kochanek & White 
(2000). 

The evolution of the collapse depends on the ratio be- 
tween the collision mean free path, A, to the local grav- 
itational scale height H (Balberg et al. 2002). Systems 
initially having ratios (A/iJ) d gc ^ 10 at the outer edge 
of the core, have a lifetime of only a few relaxation 
times 3 (Quinlan 1996; Kochanek & White 2000; Balberg 
et al. 2002). For the range of <7d m explored in Kochanek 
& White (2000), the central density increases by an or- 

1 /2 

der of magnitude in a timescale t c ,i ~ 4tT d £ n i rc , being 
Odm = 2-7rp i r s o'dm for the NFW profile. In particular, for 
Odm = 0.3, this timescale is only ~ 2.2 times the core- 
radius relaxation time. This result is troublesome for the 
survival of the core because it recollapses quickly after its 
formation and could not persist until today (see §2). 

Let us calculate the core timescale t Ci \o for the halo of 
NGC 5963. As said before, because of the adiabatic com- 
pression by the baryons, a cuspy halo, such as the NFW 
profile, is expected even in the presence of DM collisions. 
The best NFW model for T* = 1.2 corresponds to c = 15.6 
and V200 = 109 km s _1 , implying a core-radius timescale: 



3pi<r dm V 



0.8 Gyr 



Co 



1 cm 2 g _1 



(7) 



For the above halo parameters, (7d m = 2,Trpir s ad m = 
0.15 v^ 00 (ao/l cm 2 g _1 ) and hence the central density will 

(-V-ff)cdgo ~ 2 - 2 / <5 "dm- 
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be enhanced by a factor 10 in the timescale 

tcio - 4^rc = 1-25 Gyr ( J (y^) 

\100kms 1 J V lcnl S / 

(8) 

Taken a a = 0.5 cm 2 g , we get i Cj io ~ 2 Gyr. This esti- 
mate has uncertainties of as much as a factor 2. This short 
timescale indicates that the dark halo may have sufficient 
time to increase its central density by a factor 100-1000. 
In order to prevent the dark halo from such dramatic evo- 
lution, we must demand that t IC ~ one Hubble time, im- 
plying Co 0.08 cm 2 g . Proceeding in the same way we 
obtain cto < 0.05 cm 2 g -1 but now for the mass model with 
T* = 1.0. The line connecting these two points in Fig. 4 
delimits the region where cross sections are probably small 
enough to avoid catastrophic core collapse in the halo of 
NGC 5963. The corresponding line at the lcr confidence 
is also plotted. We see that the new constraint, which ac- 
counts for the role of baryons when Y* > 0.85, is tighter 
than the one derived in the previous subsections. This 
all but removes the permitted window for cross sections 
a = 0.3-1 cm 2 g" 1 in mass models with > 0.85, that 
remained open. Combining both constraints, the max- 
imum of the permitted value a ss 0.2 cm 2 s _1 , at 2er 
confidence level, occurs for Y* w 0.7. 

5. DISCUSSION AND CONCLUSIONS 

Apart from the obvious interest for the still unknown na- 
ture of dark matter, the possibility of it having a nonzero 
self-interaction cross section has other astrophysical impli- 
cations. SIDM was suggested as a route to form cores in 
LSB and dwarf galaxies. The fact that the required cross 
sections are comparable to the cross section for particles 
interacting with each other via the strong force, has led 
to speculate that DM particles could interact with both 
themselves and with baryons through the strong force 
(Wandelt et al. 2001). The interaction of dark matter 
with protons might contribute to reheat the intracluster 
medium in the central regions of clusters of galaxies (Qin 
& Wu 2001; Chuzhoy & Nusser 2004). 

A large range of parameters space is being ruled out by 
current experimental and astrophysical bounds. Previous 
studies have gradually whittled down the DM cross sec- 
tion allowed to solve the cuspy problem of halos in ACDM 
cosmology. In addition to statistical studies, the analysis 
of individual galaxies can give additional constraints on 
the strength of DM self-interaction cross section. 

The dark matter distribution in NGC 5963 is challeng- 
ing for any model designed to form central cores in dwarf 
and LSB galaxies. Apparently, its halo distribution seems 
at odds with the results of SIDM simulations with cross 
sections Odm^ioo = 0.3-1 cm 2 g _1 , which produce central 
densities ~ 0.02 M Q pc~ 3 , fairly independent of the halo 
mass, and core radii 2.5-5 kpc. The highly concentrated 
halo of NGC 5963 implies upper limits for the interaction 
of dark matter particles, unless we change our assumption 
of constant Y*. A higher Y* for the disk than for the lens 
is required to have a less concentrated halo. However, the 
observed B- V colours of the disk and the lens do not favour 
this possibility (see Bosma et al. 1988 for a discussion). 

Great efforts have been made to constrain the collisional 
cross section of DM particles. It was suggested that this 
cross section should be small enough so the core halo of 



dwarf and LSB galaxies would not collapse in a Hubble 
-l/stime (e.g., Hennawi & Ostriker 2002). Hence, the halo of 
NGC 5963 by itself cannot be going through the gravother- 
mal contraction phase. However, core collapse may be in- 
duced by the action of the baryons as they deepen the 
potential well and compress the core in an adiabatic pro- 
cess. Therefore, the core of NGC 5963 may be expanding 
if the thcrmalization has not been completed, or shrinking 
if the core is adiabatically compressed, which may occur 
only at large values of Y*. 

Either the core of NGC 5963 is expanding or contracting 
in size, a tight constraint cr < 0.2 cm 2 g _1 at 95% con- 
fidence level, which indicates that our results are robust 
to reasonable DM uncertainties, is derived. This upper 
limit for do may be overestimated as much as a factor 2 
because of our conservative assumptions. Thus, the origi- 
nal motivation for SIDM of lowering the core densities of 
galactic halos require collisional cross sections too large to 
be consistent with the halo of NGC 5963. 

One of the largest uncertainties in any mass model is 
the precise value of Y* because depends on extinction, 
star formation history, etc. Based on various considera- 
tions, e.g., stellar population synthesis, stellar counts and 
kinematics in the solar neighbourhood and kinematics of 
external galaxies, some authors argue that Y^ w 2h in 
the /-band (e.g., Mo & Mao 2000, and references therein) 
and Yf w 1.2 in the i?-band (Simon et al. 2004). This 
would mean that the constraint derived from the core col- 
lapse discussed in §4.3 is the most relevant, and suggests 
an upper limit a < 0.1 cm 2 g _1 . This value coincides 
with the upper limit inferred in clusters of galaxies (e.g., 
Meneghetti et al. 2001; Arabadjis & Bautz 2004). 

Upper limits in the 0.02-0.1 cm 2 g _1 range were derived 
by Hennawi & Ostriker (2002) from the mass of supermas- 
sive black holes in the centers of galaxies. For effective 
cross sections (7dm > 0.02 cm 2 g _1 , the accretion of SIDM 
onto seed black holes would produce excessively massive 
black holes. However, there are serious uncertainties as- 
sociated with this limit because they make use of the hy- 
pothesis that initially the density DM distribution follows 
the NFW profile, whereas the numerical study of Colin et 
al. (2002) of SIDM cosmology suggests that the NFW pro- 
file at the center of the halos is not achieved at any time. 
Therefore, the otherwise excessive accretion of matter onto 
the black hole may not occur (Colin et al. 2002). 

In order to reconcile the collisional DM hypothesis as a 
viable explanation of the formation of the constant den- 
sity cores with the halo of NGC 5963, we should identify 
potential ways able to compensate the evacuation of DM 
in the central parts caused by the scattering of DM or to 
forestall the core collapse, depending whether the core of 
NGC 5963 is expanding or collapsing. Tidal redistribu- 
tion of mass at the halo center cannot be efficient in NGC 
5963 because the nearest large galaxy in the group is at a 
projected distance of 430 kpc. One could relax our sim- 
plifying assumption that the halo consists of a well mixed 
homogeneous DM distribution and to consider a clumpy 
medium. The final density profile will be the result of 
two competing effects. On the one hand, the mass infall 
associated with the spiraling of putative massive clumps 
of DM (> 10 6 ~ 7 M Q ) towards the galactic center by dy- 
namical friction, produces a replenishment of material and 
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deepens the central potential well. On the other hand, dy- 
namical friction heating may be effective in flattening the 
inner DM profile. Unlike clusters of galaxies in which sub- 
structure is important in determining the final DM distri- 
bution (e.g., Nipoti ct al. 2004), halos of LSB galaxies are 
already largely assembled at z m 3 and hence, they have 
had sufficient time to disrupt the halo substructure, form- 
ing a smooth DM distribution within the luminous radius, 
unless NGC 5963 had a prominent mass aggregation his- 
tory. However, an anomalous accretion history would be 
very unsatisfactory because it would have dramatic con- 
sequences for the disk itself, producing an excessive dy- 
namical heating or even its destruction. Moreover, Ma & 
Boylan-Kolchin (2004) argued that energy deposition by 
merging dark matter substructures likely flatten density 
profiles. In the lack of any of those potential mechanisms 
capable to make significant mass redistribution, we think 
that our analysis is robust, even though it is based on the 
halo of a single galaxy. 

The halo of NGC 5963 is problematic for any model 
whose mechanism to produce large cores in LSB galaxies 
depends on collisions between DM particles (e.g., annihi- 
lating dark matter), reducing the parameters space and 
suggesting new directions for DM search. Our analysis 
also places strong constraints to the non-hadronic exotic 
Q-balls as a dark matter candidate in galaxies. These 
"particles" might be arranged to have mutual collisions 
with a large cross section and, in addition, Q-balls can 
stick together after collision, reducing the self-interaction 
as scattering proceeds (Kusenko & Steinhardt 2001). In 
order to be consistent with the halo of NGC 5963, self- 
interactions between Q-balls should shut off to a negligi- 
ble value after 1-2 collisions per particle in order for the 
initial scatterings to smooth out halo cusps but avoiding 



gravothermal collapse. 

Perhaps the solution of the halo core problem resides 
in decaying dark matter (Cen 200Ia,b; Sanchez-Salcedo 
2003). If dark matter particles in galactic halos decay 
to stable particles with a recoiling velocity of a few tens 
of kilometers per second, then a fraction of LSB galaxies 
can still present a substantial concentration. This novel 
dynamics associated with decaying dark matter was illus- 
trated by Sanchez-Salcedo (2003) for the case of the LSB 
galaxy NGC 3274. We believe that the halos of NGC 
5963 and NGC 3274 are not pathological cases as far as 
their DM distributions concern. There are additional in- 
dependent evidence on the existence of galaxies with large 
central densities. Loewenstein & Mushotzky (2002), in 
an as-yet-unpublishcd work, have determined the enclosed 
mass profile for the elliptical galaxy NGC 4636 using X- 
ray observations. For this galaxy the central density in 
models with dark matter cores is higher than expected in 
the SIDM scenario with a < 1 (Loewenstein & Mushotzky 
2002). 

Future observations of LSB galaxies containing a very 
low density of luminous material even in the inner parts, 
so that the observed dynamics should be dominated by 
the gravitational forces of the dark halo at small radii as 
well as large radii, will be able to determine whether cores 
are produced by the gravitational interaction between the 
luminous and dark matter, and will provide further con- 
straints on the nature of DM in galaxies. 

I thank Tony Garcia Barreto and Alberto Bolatto for 
helpful discussions. I am grateful to the anonymous ref- 
eree for constructive comments. This work was supported 
by CONACYT project 2002-C40366. 
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Table 1 

Fitting parameters for pseudoisothermal (ISO) and NFW halos after and before halo contraction. 



NFW 


ISO after decompression 


C VJoo Xred P 

(km s x ) 


r c Po 
(kpc) (M Q pc- 3 ) 


1.2 15.6 109.0 0.50 0.934 
0.7 23.2 88.0 0.44 0.962 
0.3 28.0 82.0 0.60 0.867 
0.07 33.0 76.0 0.63 0.842 


2.60 0.080 
1.20 0.260 
0.71 0.696 


ISO 




r 'c Po Xred P 
(kpc) (Mq PC" 3 ) 




1.2 0.95 0.38 0.29 0.995 
0.3 0.48 1.35 0.28 0.996 
0.07 0.42 1.75 0.28 0.996 


0.49 1.375 
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Fig. 1. — The observed CO rotation curve (empty diamonds) of NGC 5963 taken from Bolatto et al. (2003) together with the Hi rotation 
curve (empty triangles) from Bosma et al. (1988). 
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Fig. 2. — Best-fitting mass models assuming pseudoisothermal halo (left) and NFW halo (right), for different values of the stellar mass- 
to-light ratio M/L = Y* in solar units. The contribution of the luminous matter (stars plus gas) in the disk is shown as dotted lines. The 
resulting halo contribution (dashed line) to the rotation curve and the final total model curve (solid line) are also shown. The reduced x 2 of 
the fit and the probability p that the data and the model could result from the same parent distribution are given in Table 1. 
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Fig. 3. — Mass model with Y* = 1.2 assuming NFW profile for the contracted halo. Symbols as in Fig. 2. Dot-dashed line shows a fit to 
the dark matter rotation curve after adiabatic decompression (triangles). In this fit we adopted the pseudoisothcrmal profile. The parameters 
of the halo before and after adiabatic contraction are given in Table 1. 
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Fig. 4. — Exclusion contour plot for the dark matter interaction as a function of T,. Various confidence levels are shown. Solid curves 
trace the upper limit on <ro derived by the constraint discussed in §4.1. Dot-dashed lines are the core collapse (C.C.) constraint, obtained by 
requiring that the lifetime of the halo of NGC 5963 is longer than one Hubble time (§4.3). 



